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1 Introduction 



After the recent discovery of a new boson at the Large Hadron Collider (LHC) at CERN [1, 
2] one has to clarify the question whether it is the Standard Model (SM) Higgs boson or 
a new particle of an enlarged theory. For this purpose it is important to study in detail 
its production and decay processes. In the SM, but also in many extensions, the largest 
production cross section for a Higgs boson is given by the gluon fusion sub-process which, 
however, also has a large uncertainty. The uncertainties arise mainly from the parton 
distribution functions (PDFs) and from unknown higher order corrections. The latter 
amount to about 10% although perturbative next-to-next-to-leading order (NNLO) QCD 
and NLO electroweak corrections have been computed in the recent years and are thus 
available for theory predictions (see Refs. [3,4] for comprehensive reviews). In view of the 
expected experimental precision reached by the LHC experiments it is desirable to take 
the next step and compute the N 3 LO QCD corrections to the gluon fusion process. In this 
letter we contribute important building blocks to this enterprise by computing the lower 
order partonic cross sections to higher order in e = (4 — D)/2, where D is the space-time 
dimension, and by providing results for their convolution with the splitting functions. In 
our calculation we will make use of the results for the master integrals entering the NNLO 
calculation. In Refs. [5,6] results for the e expansion are provided which are sufficient for 
the N 3 LO predictions. 

The starting point for the evaluation of the cross section is the effective five-flavour La- 
grange density 

H° 

£y, e ff = qC 1 1 +C QCDi (1) 

where £qc D is the usual QCD part with five massless quarks, H° denotes the Higgs boson, 
v° its vacuum expectation value and is the matching coefficient between the full and 
the effective theory. G ^ is the gluonic field strength tensor constructed from fields and 

couplings already present in Cq CD . The superscript "0" denotes the bare quantities. 
Note that the renormalization of H°/v° is of higher order in the electromagnetic coupling 
constant. In Eq. (1) the operator 0\ is given by 

As shown in Refs. [5,7-11], through NNLO the results obtained in the effective-theory 
approach agree to good accuracy with the ones of the full-theory calculation. We assume 
that this feature also holds at N 3 LO. 

The outline of the paper is as follows: in the next two Sections we discuss in detail the 
splitting functions and the partonic cross sections up to NNLO. The procedure to compute 
the convolution integrals is described in Section 4 where also some sample results are given. 
Analytical results for all convolutions are provided in a Mathematica file which can be 
downloaded from [12]. In Section 5 we summarize our findings. 
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2 Splitting functions 



In general, the sum of the renormalized real and the virtual contributions of a given 
partonic sub-process is not finite in the limit e — > 0. The remaining poles originate from 
the collinear divergences in the initial state. For example, the initial gluon may split into 
a quark- ant iquark pair and the quark participates in the Higgs boson production. It also 
may happen that a quark in the initial state radiates a gluon before participating in the 
hard scattering process. 

The correct treatment of the singularities is achieved by the convolution of the partonic 
cross section with the splitting functions Pij(x) describing the probability of parton j to 
emit a parton % with the fraction x of its initial energy. The perturbative expansion is 
given by 

<W + ^P£'(x) +|^| 4 2) (x) + (^-) P™(x) + . . . , (3) 
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where the analytical results for P^\x) can be found in Refs. [13-16]. It is common 



practice to refer to P^ ( the fc-loop splitting function. 

We have written the expansion in Eq. (3) in terms of af \ the renormalized coupling in 
the effective theory with decoupled top quark and n% = 5 massless quarks. 

In our calculation we use the splitting functions defined in the MS scheme. This is 
consistent with the definition of the parton distribution functions (PDFs) which absorb 
the non-perturbative and non-singular features of the initial state. 

For definiteness let us present the explicit results for the one-loop splitting functions which 
read 1 
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where is the number of active flavours. As far as the two- loop splitting function Pq q is 
concerned we use P qq = + P ps [16] where the two-loop terms of and P ps are given 
in Eq. (4.6) of Ref. [15] and Eq. (4.7) of Ref. [16], respectively. The analytical result reads 
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1 Note that our definition of differs by a factor 2ni from that of Ref. [16]. 
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where ((n) is Riemann's zeta function. The remaining two-loop splitting functions can 
be found in Eqs. (4.8), (4.9) and (4.10) and the three-loop ones in Eqs. (4.14) and (4.15) 
of Ref. [16]. Note that at the three- loop level only Pgq and Pgf are needed. 



3 Partonic cross sections 

The renormalized partonic cross section at N 3 LO receives contributions from convolutions 
of splitting functions with the renormalized lower order partonic cross sections. For this 
purpose the latter needs to be expanded to higher order in e since the convolutions with 
P-j(%) come along with poles in e. 2 Furthermore, a s and the operator 0\ have to be 
renormalized which means that the LO, NLO and NNLO partonic cross sections are 
multiplied by contributions containing 1/e poles and thus the e expansion needs to be 
sufficiently deep. 

We have performed two independent calculations of the partonic cross sections. On the 
one hand we extended the full-theory calculation of Ref. [10] to one order higher in e and 
have computed the leading term in the heavy-top quark mass expansion. On the other 
hand we have implemented the effective-theory Feynman rules (cf. Eq. (1)) and evaluated 
the partonic cross sections of all sub-processes. Of course, the NLO and the NNLO results 
present in the literature [17-22] could be reproduced. The 0(e 2 ) terms at NLO and the 
0(e) terms at NNLO are new. 

From Eq. (1) it is clear that the partonic cross sections contain the two factors, the 
matching coefficient C\ and the cross section computed using the effective Higgs-gluon 
coupling contained in H O\. Thus, for the renormalized cross section we can write 

&ij = Cl&ij , (6) 

2 We refer to [17, 18] for a nice description of the procedure in case of NNLO Higgs boson production. 
In particular, in Ref. [17] one finds explicit formulae showing how the convolutions of splitting functions 
and partonic cross sections enter the renormalized result. 
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with C\ = — /(37r)(l + 0(a s )). In the remainder of the paper we will only consider a^, 
which we refer to as the "reduced partonic cross section" , since only this quantity enters 
the convolutions with the splitting functions; the multiplication with the finite factor C\ 
can be done at the very end. We note that the matching coefficient C\ has been computed 
to three- [23,24] and four-loop order [23,25,26], respectively 

In analogy to Eq. (3) we write the perturbative expansion of the reduced partonic cross 
section as 
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with A = Gfk/ (32v2) and x = M\j s. The dependence on /i 2 /M|, where /i represents the 
renormalization and the factorization scale, is suppressed. At NLO we have i,j G {g, q, q} 
and at the NNLO the possible initial states are given by gg, qg, qg, qq, qq, or qq', where 
q and q 1 stand for (different) massless quark flavours. 3 The computation of the N 3 LO 
contributions requires that is known to order e 3 ~ fc . 

In order to fix the notation we provide the explicit results for af^ and a^j which are given 
by 
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(8) 



with x^h — ^/Mjj. The other results for the partonic cross sections are significantly 
larger and are given in computer readable form in Ref. [12]. 



4 Convolutions of partonic cross sections and split- 
ting functions 

The convolution of functions / and g that enter P^\x) and bij is defined as 



[/ <g> g] (x) 



dxidx 2 5(x - xix 2 )f(x 1 )g(x 2 ) 



(9) 



3 It is understood that ghosts are always considered together with gluons. 
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Both the results for the partonic cross sections (see previous Section) and the splitting 
functions include combinations of HPLs up to weight four with factors 1/x, 1/(1 — x), 

and 1/(1 + x), and the generalized functions 5(1 — x) and ln ^j— — 

briefly describe our approach to compute the convolutions of such functions. A detailed 
description can also be found in Appendix B of Ref . [5] . 



. In the following we 



1. First, we transform the convoluted expressions to Mellin space via 



M n [f(x)} 



dxx n -'f(x). 



This transformation turns convolutions into products due to the identity 



M„ 



[f ® g] (x) 



M n [f(x)]M n [g(x)} . 



(10) 
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A comprehensive discussion of the Mellin transform and the list of all Mellin images 
appearing in the calculation of the NNLO Higgs boson production rate can be found 
in Refs. [27,28]. Mellin transforms of HPLs and their derivatives can be conveniently 
expressed in terms of harmonic sums [29,30]. In our calculation, we used the FORM 
package harmpol [31] to perform the transformation. 

2. Second, we prepare a table of Mellin transforms of HPLs to certain maximum weight. 
For example, the table through weight one reads 



M n [l] 
M„[H (x)] 
M n [H!(z)] 
M n [H_!(x)] 



1 

n 



Si(ra) 



n 



-if/a / \ , oN In 2 

n n 



(12) 



where Hj(x) are the HPLs of weight one and Sj(n) are the harmonic sums of weight 
one. Mellin transforms of HPLs of higher weights produce harmonic sums of higher 
weights and transcendental numbers originating from HPLs evaluated at x = 1. 

3. Third, we prepare a table of Mellin transforms of "regularized derivatives" of HPLs 
in the table above as follows: 



M„ 



M„ 



OA 



d x f(x) =R[f(x)]-(n-l)M n _ 1 [f(x)} 



where 
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(13) 
(14) 

(15) 
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with gj (l) ^ Vj > 0. 



This definition allows us to establish relations between the regularized derivatives 
of HPLs and the "common" generalized functions. In particular, 
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Note that if f(x) is not divergent at x — 1, its regularized derivative Eq. (14) reduces 
to the usual derivative: 



ax 



x n - 1 f(x) Q -(n-l)M n ^[f(x)]. 



(17) 



For example, 



4H_i(x) = ^-H_!(x) = -J— , 

cix L -\- x (18) 

4,H_i,i(x) = ^-H_ 1:1 (x) = 

ax 1 + x 

Thus, it is not necessary to separately consider Mellin transforms of all functions 



appearing in the derivatives of HPLs, such as 
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and so on. Rather it 



is sufficient to consider Mellin transforms of regularized derivatives in case Mellin 
transforms for generalized functions are needed. For illustration we list the Mellin 
transforms of regularized derivatives of HPLs to weight one 

M n [d x l] = 1, 
M n [4H (x)] 
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4. Finally, we prepare and solve the system of linear equations and thus perform the 
inverse Mellin transform. Combining the tables for Mellin transforms of HPLs as 
in Eqs. (12) and the regularized derivatives of HPLs as in Eqs. (19) one obtains a 
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system of linear equations for terms of the form l/n k , S...(n)/n fc and (—l) n S...(n)/n k 
which can be used in order to find the inverse Mellin transforms. In principle, the 
expression in the Mellin space may have other terms, whose inverse transform cannot 
be determined from the system, but all such terms cancel in the real calculation. 

The above steps have been implemented in a Mathematica program which uses the HPL 
package [32, 33] and can compute all the necessary convolutions. 

As an important cross check, we have performed all the convolution integrals also numer- 
ically for some specific value of x. In case no plus distribution is present in the integrand 
it is straightforward to perform the numerical integration. For convolution integrals in- 
volving one or more plus distributions, the cancellation of the singular behaviour among 
terms can be quite involved. The problem can be solved by introducing auxiliary regu- 
larizations [34] 

1 } lim lim —^-(1 - xY l+a \ (20) 
1- x v->o 7] n da n 

for evaluating integrals of the individual terms. After performing the integrals by using 
S functions in the definition of the convolutions one rescales the integration variables 
such that the integration domain becomes the unit hypercube. The singularities come 
from both the lower and the upper limit of the integration as well as from the overall 
(1 — x)~ 1+ari factors which give 5 functions and plus distributions in the result of the 
convolution, after applying 



arj — J! 



ln j (A) 
A 



(21) 



In the next step one differentiates the result with respect to a, takes the limit a — > 1, and 
finally expands in 77. All pole terms l/r] k should cancel among the integrals and hence 
the leading term gives the result of the convolution. The remaining integrals have no 
divergences and can be performed numerically. 

As a further welcome check we have compared our results with the explicit expressions 
for the convolutions given in the Appendix of Ref. [35] . In particular we found complete 
agreement with Eqs. (C.28) to (C.31) which involve convolutions of two plus distributions. 

In Tab. 1 we classify the convolutions which are required for the N 3 LO calculation accord- 
ing to the number of nested convolutions. Contributions where only one convolution is 
involved already enter the NLO and the NNLO results. At the N 3 LO, the LO, NLO and 
NNLO partonic cross section has to be convoluted with three-, two- and one-loop splitting 
functions, respectively. Two convolution integrals are only present at NNLO and N 3 LO 
and three convolutions only at N 3 LO. In the latter case, however, only one-loop splitting 
functions and the leading order cross section enter. Since the latter is proportional to 
5(1 — x), only two non-trivial integrals remain. 

The right-most column in Tab. 1 defines the name of the file where the analytical results 
can be found. At NLO there are only 2 non-zero convolutions, at NNLO there are already 
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convolution 


order 


file name 
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StilO<ij>-PKjk>-Pl<kl>-Pl<lm>.m 
Pl<im>-StilO<ij>-PKjk>-PKkl>.m 



Table 1: Classifications of the convolutions needed up to order N 3 LO where k,l,m G 
{<?> 9) 9) ?'}■ As far as splitting functions are concerned there is no difference between "g", 
"g" and u q' n . Furthermore, we have = <7-™\ 



18 and at N 3 LO altogether 82 convolutions of partonic cross sections with up to three 
splitting functions. All the results can be downloaded from [12] where for convenience we 
also provide intermediate results for the convolutions involving only splitting functions. 
In Appendix A we provide a list assigning the convolutions to the file names. Note that 
not all results are shown which are identical due to symmetry. 

The graphical representations of the five different types of convolutions given in Tab. 1 
are displayed in Fig. 1. 

Let us at the end list some convolutions which appear as building blocks in the course of 
our calculation. 
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Figure 1: Schematic representation of the convolutions of the partonic cross section 5"^- 
with the splitting functions Pjk- 

+ yC(3) - 3C(5) . (22) 

5 Summary 

In this paper we computed the NNLO partonic cross sections for the Higgs boson produc- 
tion to order e which constitutes a building block for a N 3 LO calculation. Furthermore 
an algorithm has been developed which allows the evaluation of all convolutions of a 
partonic cross section with splitting functions needed at N 3 LO. The described proce- 
dure might also be interesting for other processes like Drell-Yan production. We provide 
computer-readable results [12] both for the partonic cross sections and for all the con- 
volutions entering the N 3 LO calculation of the partonic cross section for Higgs boson 
production at the LHC. 
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A Convolutions and names of result files 

Tables 2 and 3 list the convolutions entering our calculation together with the names of 
the result files [12]. For convenience we use the notation q = b and q' = p. 
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Table 2: Correlation of convolutions involving the partonic cross section and names of 
result file. 
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Table 3: Correlation of convolutions involving only splitting functions and names of result 
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